CU-TP-776 



Dependence of QCD hadron masses on the number of dynamical quarks 

Dong Chen^* and Robert D. Mawhinney^* 

'^Department of Physics, Columbia University, New York, NY 10027, USA 

We have studied the hadron spectrum while varying the number of light dynamical quarks when the physical 
lattice spacing and volume are held fixed relative to the rho mass. For two and zero fiavors of staggered fermions, 
we find the nucleon to rho mass ratios (extrapolated to zero valence quark mass) are very similar. However, for 
four flavors the ratio is 7% (2 a) above the two flavor result. 



1. INTRODUCTION 

The effects of dynamical quarks in zero tem- 
perature simulations of QCD have generally been 
found to be quite small for ratios of physical quan- 
tities. In particular, quenched and two flavor dy- 
namical simulations, for similar lattice spacings 
and volumes, have yielded similar results for the 
nucleon to rho mass ratio. Currently we can not 
probe the differences between zero and two flavor 
calculations as the continuum limit is approached, 
due to limitations in computer power. Instead, we 
have undertaken a comparison of quenched, two 
flavor and four flavor QCD at zero temperature, 
while holding the lattice spacing and volume fixed 
in physical units, to look for dynamical quark ef- 
fects. 

At Lattice '95, the Columbia group reported 
a spectrum study with two flavors of staggered 
quarks on a 16"^ x 40 lattice 0|. Since then, sim- 
ulations for four and zero flavors on 16'^ x 32 vol- 
umes at Columbia and zero flavors at Ohio State 
have been done. We have chosen our simulation 
parameters so that the valence rho mass, extrap- 
olated to zero valence quark mass (rrivai = 0), is 
independent of the number of dynamical quarks. 
(We have achieved this to the 3% level.) All the 
simulations consist of large data samples (by cur- 
rent standards), affording a precise comparison. 

*The 2 flavor calculation in tliis work was done in collabo- 
ration with Shailesh Chandrasekharan, Norman H. Christ, 
Weonjong Lee, and Decai Zhu; The 4 and flavor calcu- 
lations were done in collaboration with Norman H. Christ 
and Gregory W. Kilcup at the Ohio State University. This 
work was supported in part by the Department of Energy. 
Presented at Lattice '96. 



2. SIMULATION DETAILS 

Table |l| lists the parameters for our simula- 
tions. The Columbia simulations were all done on 
the 256-node, 16 Gigaflop computer at Columbia, 
now completing its seventh year of full-time calcu- 
lations. The two flavor simulation was done with 
the inexact R hybrid molecular dynamics algo- 
rithm of ; the four and flavor simulations used 
an exact hybrid Monte Carlo algorithm, employ- 
ing the $ hybrid molecular dynamics algorithm 
of [|| and a Monte Carlo accept/reject step. The 
quenched simulations of OSU were done on the 
Ohio Supercomputer Center T3D using a mixture 
of over-relaxed and Metropolis steps. 

From Table |^ it is apparent that the two and 
four flavor simulations had very similar param- 
eters for the evolution. For the four flavor cal- 
culation, we used a tighter conjugate gradient 
stopping condition in the evolution, so that the 
accept/reject step of the exact algorithm would 
be based upon as accurate a value for the 5- 
dimensional Hamiltonian as possible. The stop- 
ping condition we used, and a test condition a 
factor of 3 smaller, both gave the same sequence 
of lattices for about 10 time units when evolved 
from the same starting lattice. 

The four flavor run had an acceptance rate of 
95%. This high acceptance rate is reassuring for 
the choice of the running conditions for our inex- 
act, two flavor run. In particular this is evidence 
that the errors in our two flavor run due to using 
an inexact algorithm are small. This is impor- 
tant, since we want to compare physics from the 
two different algorithms. 
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Table 1 

Simulation parameters for the four runs presented. The run length, thermalization, hadron measurement 
frequency and jackknife block size are in time units for the CU runs and sweeps for the OSU runs. 





Nf = 4 (CU) 


Nf^2 (CU) 


Nf = (CU) 


Nf^O (OSU) 


V yji. uii-Livj 


16^ X 32 


16^ X 40 


16^ X 32 


16^ X 32 




5.4 


5.7 


6.05 


6.05 


^dynamical^ 


0.01 


0.01 






evolution 


HMC 


HMD 


HMC 


OR + Metropohs 


run length 


4450 


4870 


187,125 


787,500 


thermalization 


250 


250 


375 


25,000 


acceptance rate 


0.95 




0.91 




trajectory length 


0.5 


0.5 


0.75 




step size 


0.0078125 


0.0078125 


0.025 




CG stopping condition 


1.13 X 10"" 


1.01 X 10"" 






total run time 


5 months 


7.5 months 


1.7 months 


25k node-hours 


hadron source 


16"^ wall, all coordinates {x,y,z) even 


valence quark masses 


0.01, 0.015, 0.02, 0.025 


hadron measurement frequency 


5 


6 


187.5 


2500 


number of lattices 


840 


770 


996 


306 


measurements per lattice 


4 


5 


1 


4 


jackknife block size 


50 


60 


2250 


7500 


number of blocks 


84 


77 


83 


102 



3. RESULTS 

Figure |^ shows pi, rho and nucleon effective 
masses for the four simulations, where the effec- 
tive masses are found from fits to time slices t to 
t + 3. The wall source is fixed to Coulomb gauge 
and the simple local staggered sinks are used. 
While the rho masses for all cases are very similar, 
the four flavor nucleon mass is significantly larger 
than the others. Notice there is little difference 
between zero and two flavors. 

Figure || shows the dependence of the rho and 
nucleon masses (found from fits between t = 6 
and t = Nt/2) on m^^\. The rho masses agree 
quite well for all values of rnvai, not just for 
TOvai ^ 0. A consistently larger value for the nu- 
cleon mass for four flavors is apparent. We have 
tried breaking our two flavor data into halves and 
re-fitting the masses for each half. The resulting 
errors are consistent with the expectation that 
they be larger. Similarly for four flavors we 
compared results from the first half of our sample 



with the full sample and the errors are consistent. 
Since we see no evidence for grossly underestimat- 
ing our errors for half of our samples, we believe 
the errors shown in Figure |^ are reasonably reli- 
able. 

Figure ^ shows the dependence of the nucleon 
to rho mass ratio on Wvai- The lines are ex- 
trapolations of the ratios, while the data points 
for TOvai = are the ratio of the extrapolated 
masses. The ratio of the extrapolated masses 
gives TOjv/mp equal to: 1.591(41) for Nf — 4, 
1.489(31) for Nf = 2, 1.470(29) for CU iV/ = 
and 1.443(28) for OSU Nf = 0. Note the 1 cr 
agreement between the two quenched calculations 
and the closeness of the two flavor result. The 
four flavor result is 2 cr above the two flavor re- 
sult. 

The error on the four flavor ratio of mi^lm,p 
is larger than for two flavors, while the error on 
the individual masses is quite similar. We have 
studied this in some detail and observe that for 
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Figure 1. The effective masses for the tt, p and A'' 
for mvai = 0.01. 



four flavors, for the block size we are using, there 
is less correlation between fluctuations in the rho 
and nucleon propagators than for the other simu- 
lations. When the ratio is calculated, there is less 
cancellation between these uncorrelated fluctua- 
tions, leading to the larger error quoted. 

In conclusion, we have seen a 7% difference (2 
a) in viM /iTip for rrivai = depending on whether 
two flavors or four flavors of dynamical fermions 
were used in the evolution of the configurations. 
Notice that the increased number of flavors is 
making the unphysically large value for the ra- 
tio increase. We are anxious to understand if this 
effect persists for larger volumes and when the 
four flavor dynamical fermion mass is varied. 

We are currently checking for the effects of four 
dynamical fermions in other observables. 



Figure 2. nip and mjv extrapolated to mvai 
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Figure 3. mN/mp versus rrivai- The lines are 
extrapolations in the ratios and the data points 
at mvai = are ratios of extrapolated masses. 



